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The coupling coefficients of the modes of a parabolic index fiber with 
randomly curved axis are derived and are used to compute its excess losses 
and impulse response. It is found that bends with a period comparable 
to the natural ray oscillation period in the parabolic index medium are 
catastrophic. The average radius of curvature R c of a guide composed of 
circular sections with an average length of 1 cm must not decrease below 
approximately R c = 1 m. Mode coupling by random bends has the tend- 
ency to reduce the width of the impulse response function. However, this 
improvement is accompanied by losses. Reducing the width of the impulse 
response for coupled mode operation to half its uncoupled width causes 
0.7 dB additional loss, a ten-fold reduction of the pulse width costs 18 dB. 

I. INTRODUCTION 

Optical fibers with parabolic refractive index profile 1 ' 2 have less 
pulse delay distortion than conventional fibers with piecewise con- 
stant, discontinuous index distribution. 3 The width of the impulse 
response increases in direct proportion to the length of the fiber. It 
is well known that an improvement of the impulse response results if 
the modes are coupled among each other. 4 - 5 In the presence of mode 
coupling the width of the impulse response increases only propor- 
tionally to the square root of the length of the fiber. 

Pulse propagation in multimode parabolic index fibers is studied in 
this paper by means of converting the coupled power equations to a 
partial differential equation. 5-7 Random changes of the direction of 
the waveguide axis are considered as the coupling mechanism. 

The problem is simplified by assuming that the modes of the para- 
bolic index fiber are essentially the same as the modes of an infinitely 
extended square-law medium. The fiber boundary is included in the 
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description by requiring that modes interacting with it suffer high 
losses so that an effective cutoff exists. Modes below the cutoff value 
propagate as if they were in an infinitely extended medium. At cutoff we 
demand that the modes do not carry power. The effect of the wave- 
guide wall is thus taken into account as a boundary condition that 
has to be satisfied by the solutions of the partial differential equation. 
This formalism provides information about the width of the impulse 
response function and the losses associated with the coupling mecha- 
nism. The achievable improvement in pulse width due to mode coupling 
can thus be expressed in terms of the associated loss penalty. We 
conclude that mode coupling is capable of improving the already 
favorable impulse response of the parabolic index fiber. However, 
this improvement of the width of the impulse response is accompanied 
by excess losses. The product of the square of the pulse width ratio 
(width of the impulse response of coupled modes to the uncoupled 
pulse width) times the loss penalty is independent of the waveguide 
parameters and the statistics of the axis deformation. There is thus 
no hope of reducing the loss penalty of delay distortion improvement 
by optimizing the waveguide parameters. 

II. MODES AND COUPLING COEFFICIENTS 

We use the modes of the infinite square-law medium. The refractive 
index distribution is assumed to be of the form 

n' = ng(l-2A^)- (D 

The fiber radius is at r = a. However, the modes are assumed to be 
unaffected by the fiber boundary if their mode number remains below 
a cutoff value. We use linearly polarized modes and obtain for the y 
component of the electric field 8 

2 (+fep)* H v (^2-\H q (rt V -\ <r r2 "" 2 

Mjpq (n 7r2p+«p!g!)ny 

There is also an electric field component in axial direction. But, for 
small refractive index changes, it is negligible. The functions H p and 
H q are Hermite polynomials of order p and q, the radius r is defined by 

r 2 = X 2 + y 2 } (3) 

and the mode radius w is 

V2a V (4 ) 



w 
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with the free-space propagation constant 

k = oWeoMo. (5) 

The propagation constant of the mode is given by the expression 8 

*- B « k [ i -S (p+s+i) T- (6) 

The orthogonality of the modes and their normalization follows from 
the following equation : 



2k 



- f r E pq E' p . g .dxdy = P8 pp .8 gq .. (7) 

Mo J— to J—x, 



The asterisk indicates complex conjugation. 

The cutoff condition for the guided mode has been derived in Ref. 9. 
The permitted maximum values of p and q are defined by the relation 

[K a 2 

(v + q)* = Jg n ° ka = w 2 ' ^ 

Any deviation of the parabolic index fiber from its perfect geometry 
can be expressed by a change of its refractive index distribution. 

Changes in the direction of the waveguide axis can be expressed by 
the following index distribution : 



n 2 = nl 



i -2j[(*-/d))» + i«] 

I* 



(9) 



We consider waveguides bends in only one plane for simplicity. The 
results thus obtained can easily be extended to the general case. The 
appropriate coupling coefficient for this type of index change is 1<U1 

Kpq - p,Q ' = 4P(0,r%,v) /-" /-" ^ ****'«•**■ ao) 

Bends of the waveguide axis couple even modes to their immediate 
odd neighbors and odd modes to their immediate even neighbors. 
Only the following coupling coefficients are different from zero : 

df 

v n kwA r dz . n kwA i— , . , lt . 

K P9lP± i, 9 = —-T- Vp = i —^— ipf(z). (U) 

a P P Pp±\ a 

Because we restricted the waveguide curvature to the x — z plane 
there is no coupling between the modes with different values of q. 
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All coupling coefficients with different q values vanish. The derivative 
of f(z) was replaced by the function itself with the help of the relation 

1 »¥-></«. < 12 > 



- 0' dz 

This replacement is permissible since it is the spatial frequency ft — p' 
of f(z) that is responsible for the coupling process. 

III. COUPLED POWER EQUATIONS 

Pulso propagation in multimode optical fibers can be described by 
the following set of coupled equations for the average power P, 
carried by the modes : 5 

tt* + ~ ^r - - a * p * + £ *"< p " - p " } - (13) 

dz Vf, at v = i 

The single label v indicates both p and q. The power coupling coeffici- 
ents h vii are defined by 5 '° 

K, = h,„ = \&r\*F(fi,-fi,). (14) 

The coupling coefficient (11) enters the power coupling coefficients 
via the definition 

Kp, = £„/(*). (15) 

The power spectrum F (£„ - /?„) is defined by the equation 

F{6) ={14= / o L /W«-*a|\ (16) 

The symbol ( ) indicates an ensemble average. 

Coupling between the modes of the parabolic index fiber is an ideal 
application for a diffusion theory of the power coupling process. 7 
Since only nearest neighbors couple directly to each other, power 
redistributes itself by jumping from mode to mode in the same way 
as particles diffuse through real space. If the mode number is very 
large we can consider the set of discrete modes as a quasi-continuum 
and change the equation system (13) into a partial differential equation. 
To accomplish this transformation we consider the following expres- 
sion using h^y = h, 






Z K,(P, - P„) = Am+lmCPm+i - P,) - K^iP, - P„_i). (17) 
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Considering y. as a continuous variable y. = 6 we use the approximation 
A(i+i,x(P/i+i — i%) — hp^-iiP,, — P M -i) 

-H" + -»(SL.-™(5).1 

With Ad = 1 we can write (13) as a partial differential equation 

= - p + £(*£)■ (19) 



a? lap 

d2 v dt 



The propagation constant (6) can be approximated as follows : 

= n k - — (0 + q + 1) (20) 

Lv 

with p = 6. The difference of the propagation constants of adjacent 
modes, 

A/3 = 0(0 + A0) - 0(0) = - — = - Q, (21) 

is independent of 0. The power spectrum entering (14) contributes to 
the coupling process only at one fixed spatial frequency and is in- 
dependent of the variable 0. The power coupling coefficient (14) can 
be expressed with the help of (4), (11), and (21) as follows (p = 0) : 

A W .*5*l , m ,. (22) 

With (22) the partial differential equation (19) assumes the form 
dP , 1 dP n , r [J J P , afl 

Tz+vl>i=- aP + K [ d W + ire\ (23) 

with 

K = £g V * < fl >- (24) 

a 3 

We assume that the attenuation coefficient a in (19) is constant and 
describe the high loss, that we must attribute to modes interacting 
with the waveguide boundary, by means of the boundary condition 

P(z, t, e) m o for = e e . (25) 
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The cutoff value 6 = 6 C follows from (8) : 

IY 

- q. (26) 



a 2 _ WokavA 

'-' v ~ q vT 



vr 



The slope dP/dd determines the rate of power diffusion. Since no 
power can be lost at = we must also require 

^ = at = 0. (27) 

do 

IV. STEADY-STATE POWER LOSS 

We begin the discussion of the solutions of (23) by neglecting the 
fact that each guided mode has a slightly different group velocity and 
consider v(B) = const. We construct a solution of (23) by introducing 
the trial solution 

P(z, t, 6) - e-^ +a) 'G{e). (28) 

Substitution of (28) into (23) yields the ordinary differential equation 

•SF+£ + §°- - (29) 

The normalized solutions of this equation that satisfy the boundary 
conditions (25) and (27) are 



G,(0) = -i 



! /.( 



U»\hr 



(30) 



4e c JiM 

with 



w > =2 V^- (3i) 

The parameters u v are determined as the roots of the equation 

/,(«,) = 0. (32) 

The functions G v {6) are mutually orthogonal. 

The general solution of the power equation (23) is obtained as the 
superposition of the trial solutions 

P(z, t, d) = e-°> £ cGMer'". (34) 

1=1 

The expansion coefficient c, can be determined from the power dis- 
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tribution at z = with the help of the orthogonality condition (33), 

c, = I ' G,(0)P(O, t, 6)dd. (35) 

Jo 

The eigenvalues a v are obtained from (31) and (32), 

°'-W e ' (36) 

The eigenvalues increase with the increasing values of the roots w„. 
It is thus apparent that only the first term in the series (34) needs to 
be considered for large values of z. The steady-state power distribu- 
tion is thus described by the equation 

P(z, t, 6) = Cie- (a+ffl K?i(0) for z -> oo. (37) 

After an initial transient has decayed, the power distribution (versus 
mode number 6) assumes the steady state (37). The power loss in the 
steady state is the sum of the constant loss a, that was assumed to be 
the same for every mode, plus the loss value <ti that stems from mode 
coupling due to waveguide curvature. With 

«i = 2.405 (38) 

we obtain the steady-state curvature loss from (24), (26), and (36), 

_ 2.045 -njtg! 

" " T^m r, my (39) 

Because of our assumption that the waveguide is curved only in one 
plane the steady state losses depend on the mode number q. For small 
values of q we find low curvature loss 

T»=-^f(a) for q = 0. (40) 

With increasing values of q the losses increase until they reach in- 
finitely high values. 

In any actual cases it is unrealistic to assume that the waveguide 
would be bent in only one plane. Bends in the perpendicular plane 
couple modes with different q values. The total steady-state loss is 
thus a weighted average of the losses (39). The weight factor is the 
number of modes for each value of q. According to (8) we have 

_ /I 



p = N(q) = Jtlnoka - q (41) 
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different modes for each value of q. The average loss that results from 
coupling all the modes by random bends in both planes is thus (see 
appendix) 



o\ = 



(noka) 



/•nokaVA/2 5 8 A , , 

* f <nN(q)dq = ii -^F(a). (42) 

Comparison with the loss coefficient (40) for the mode group with 
q = shows that the total fiber loss (42) is just twice as large. The 
loss coefficient (40) is representative of a slab waveguide. The actual 
loss of the round fiber can thus be deduced from the slab waveguide 
model. That the loss coefficient of the fiber is twice as large as the slab 
waveguide loss might be expected, since F(Q) stands for the ampli- 
tude of the power spectrum for bends in only one plane. The fiber is 
assumed to be bent in both planes with equal power spectra. The 
effect of both bends add, doubling the loss coefficient. 

In a fiber cable the fibers may (or may not) be twisted around each 
other. Such twists could introduce an almost sinusoidal deformation of 
the fiber axis. The length A of the period of sinusoidal deformations 
that should be avoided follows from (21), A = 27r/fi. For numerical 
estimates we are using a fiber with radius a = 4.85 X 10" 3 cm and 
A = 0.014. The critical period for this fiber is A = 0.18 cm. If such a 
period should be built into the fiber by the method of cable construc- 
tion we can estimate the losses that would be caused by a given ampli- 
tude. F(S2) has the dimension of cm 3 . It can be interpreted as the ratio 
of the square of the amplitude of the sinusoidal deformation and the 
spatial bandwidth that may be caused by random phase changes. 
With U = 34.5 cm -1 let us assume a spatial bandwidth of AQ = 3 
cm" 1 . An excess loss of a x = 10 dB/km = 2.3 X lO" 5 cm" 1 would re- 
quire F(fl) = 1.57 X 10- 13 cm 3 . The square of the amplitude A is 
given by the product of F (12) with the spatial bandwidth. We thus 
obtain the amplitude of the sinusoidal deformation of the fiber axis 
that causes an excess loss of 10 dB/km: A = [F(a)AO]* = 6.9 
X 10 -7 cm = 69 A. Sinusoidal axis deformations at the critical wave- 
length are thus seen to be extremely dangerous. 

V. LOSSES FOR A STATISTICAL MODEL 

Even though it is only one spatial frequency of the power spectrum 
F(0t — (SO that determines the steady state loss (42), it is hard to 
guess the amplitude F(Q) that might be expected at this spatial fre- 
quency fi. In order to gain insight into the expected steady-state 
curvature losses it is necessary to consider statistical models. We 
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consider a model consisting of waveguide sections with constant 
curvature whose magnitude and sign varies randomly. 
The power spectrum can be written as follows : 

l /I r L l 2 \ 

f(q) = i\\j Q /«*-"■*!> 

VL\\J, dz* e I/- (4dj 

The step from the function f(z) to its second derivative involved two 
partial integrations. The end points of the integration range do not 
contribute if we assume that the randomly disturbed guide is con- 
nected to two perfectly straight waveguide sections so that f(z) and its 
first two derivatives vanish at z = and at z = L. 

For waveguides that are only slightly bent we can consider the second 
derivative of f(z) as the curvature function 1/R c (z). We denote with 
C(u) the autocorrelation function of the curvature function, 

0(u) = ( p,*, , — s V (44) 

\R c (z)R c (z + u) / 

It is well known that the power spectrum of a function is equal to the 
Fourier transform of its autocorrelation function. 12 We may thus write 

F{$) - I £ C{u)e-**«du. (45) 

The autocorrelation function of waveguide sections with piecewise 
constant curvature and fixed length D is 

C(u)=\ D K \ U \= L > (46) 

[o \u\ > D. 

The parameter k 2 is the variance (square of the rms value) of the curva- 
ture l/R c . Substitution of (46) into (45) results in 

F = m« (1 " cos QD) - (47) 

An average over this expression, that allows us to consider D as an 
averaged quantity, leaves us with 
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From (21), (42), and (48) we obtain the following expression for the 
steady-state loss of our statistical waveguide model : 

*<■> = *&. m 

As a numerical example, we consider a parabolic index waveguide 
with the following parameters : 

X = 1 /zm, free-space wavelength (k = 6.28 X 10 4 cm -1 ) ] 
a = 4.85 X 10~ 3 cm, waveguide radius 
n = 1.56 
5 = 0.014 
w = 7.69 X 10 -4 cm, mode radius 

^ = 6.3 

w 

With these data we have (k in cm -1 , D in cm) 



(50) 



<« = 0.17 jj cm- 1 . (51) 



We may now ask for the rms value of the curvature that is required 
to cause a steady-state loss of 2.3 X lO" 8 cm" 1 = 10 dB/km with an 
average length of the waveguide sections of D = 1 cm. We find from 
(51) k = 0.0116 cm" 1 or 1/k = 86 cm. Our result tells us that a wave- 
guide composed of individual sections of constant curvature of average 
length D = 1 cm with R c « 1 m radius of curvature has 10 dB/km 
additional loss. For our derivation we assumed that the high-order 
modes suffer very high losses since their fields reach into the vicinity 
of the waveguide wall. Whether the interaction with the outer wave- 
guide boundary causes high losses depends on the construction of the 
waveguide. If the outer surface is rough or coated with an absorbing 
material to reduce crosstalk, the losses are high and our estimate 
applies. 

VI. PULSE DELAY DISTORTION 

It was shown in Ref. 5 that the width of the impulse response of a 
multimode fiber with coupled modes is given by the equation 

M = 4VpL. (52) 



L is the length of the waveguide and p is the second-order perturbation 
of the eigenvalue <n defined by (36). For the discrete case we write 
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G y (6) = G p <- V) and obtain p in the form 5 



N 



LMdh'H'. 



.=2 tr, — ffi 



The average group velocity v actually does not contribute to (53) 
on account of the orthogonality of the vectors G {>,) . With the assump- 
tion of a continuum of modes we obtain instead of (53) 



p= z 



^U!\w)-i) 6im ' im J\' (54) 

The inverse group velocity is obtained by approximating (6), 

0« n k - ^? (p + q) - -jL (p + 9 )2, (55) 

and taking the derivative. With v = c/n we obtain (with p = 8 and 
c = light velocity in vacuum) 

mri^£-*-?B9<'+»i+& (56) 

The term with q 2 does not contribute to the following integral because 
of the orthogonality relation (33) : 

- 16 gw ,"-"v.<4 i -T^ ± ^1+»l- ^ 

cn k 2 a 2 (m; — wf) 2 [ (_ (w 2 — w? ) 2 J ] 

Each mode group with a given value of q has a different spread of 
the group velocities of its uncoupled modes. However, we have seen 
that the waveguide losses could be obtained from the simpler slab 
waveguide model. This simplification is expected to apply also to the 
pulse distortion problem. The slab model is obtained by setting q = 0. 
The spread of inverse group velocities is largest for q = since the 
allowed range is largest in this case. However, even though the spread 
of the group velocities is reduced for increasing values of q, the mode 
groups with different q values arrive at different times. This delay 
distortion is reduced by coupling of the different mode groups by 
means of waveguide bends in the perpendicular plane (perpendicular 
to the plane coupling the modes with different p values). The mode 
group with q = can be excited by shining light into the fiber that 
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is collimated in one plane but spreads in the plane of the bends in 
such a way that all modes with 5 = and p values in the range 
< p < n ka(A/2)* are excited. The bends in one plane do not 
cause coupling to modes with different q values but reduce the delay 
distortion of the modes with different p values. From this physical 
picture we see that the delay distortion problem is reduced to studying 
delay distortion in a slab waveguide. Bends of the fiber in the perpen- 
dicular plane couple the modes with different q but fixed p values. 
Their velocity spread is the same as that considered in the first problem. 
We thus expect to obtain the correct result by considering the delay 
distortion reduction for the mode group with q = 0. 
With q = we obtain from (26) and (57) 



C{w)-i) Gl{e)G ' (e)de 



8n„A 2 u v Ui 



t-m + T&l (58) 



c (ul - u\Y _ (ul -u\Y 

Using (24), (36), (54), and (58) we have 

12871*0,* A 3 ( ,=. u\u\ [\ \2(ul + u\) 



P = 



- u\u\ [ _ 

hiui-uirl 1 - 



c*F(Sl) I ,t2 (ul - w?) 5 L M ~ "?) 2 . 



= 2.26X10-.^A- (59) 
The width of the impulse response follows from (52), 

From Ref . 9 we obtain the width of the impulse response for uncoupled 
modes 

At = ?£ 5». (61) 

2c 

The improvement that is caused by mode coupling is the ratio of the 
widths of these two impulse responses 

D M °- 12 ° 2 (RO\ 

R = j- T = m®m ' (62) 

For the statistical model of a sequence of circularly bent waveguide 
sections we obtain with (48) and (21) 

0.24A /D\* 



R - 



(*)'• (63) 
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We may ask for the average radius of curvature that is required to 
cause a ten-fold improvement of the width of the impulse response 
due to mode coupling. Using L = 1 km and an average length of the 
bent sections of D = 1 cm and the numbers in (50) we obtain for 
R = 0.1, k — 0.022 cm -1 or an average radius of curvature R = 1/k 
= 45 cm. This relatively small radius of curvature may cause very 
substantial excess loss according to the loss example in Section V. 

In order to relate the excess loss to the delay distortion improvement 
we consider the loss penalty that must be paid for a given amount of 
improvement in the width of the impulse response. Both the loss 
formula (40) and the improvement factor R, (62), contain the power 
spectrum of the distortion function /(z). By taking the square of R 
and multiplying it with the loss per length L, o\L, we obtain t 

R 2 a iL = 0.042 = 0.18 dB. (64) 

This important formula is independent of any of the waveguide 
parameters and of the statistics of the axis deformation. This means 
that the loss penalty, a\L, for parabolic index fibers depends only on 
the delay distortion improvement that one wants to achieve. For 
R = 1 we have a loss of o\L = 0.18 dB. Clearly, the range of ap- 
plicability of (64) is exceeded in this case since R = 1 means that 
there is no improvement at all. For R = 0.5 we pay a loss penalty of 
0.7 dB, 22 = 0.1 increases the loss to 18 dB. The already favorable 
delay distortion of the parabolic index fiber can be improved by 
intentional curvature of the waveguide axis. 

VII. DISCUSSION 

We have studied the performance of the parabolic index fiber with 
randomly curved axis. The curvature of the waveguide axis has the 
tendency to force a light beam inside of the fiber towards the fiber 
boundary. In terms of wave optics this means that the wave field 
begins to interact with the boundary of the fiber. If this boundary 
is perfectly smooth no particular harm may be done except that the 
impulse response of the fiber is likely to deteriorate. However, the 
interfaces between two dielectric regions tend to be rough. Surface 
roughness leads to scattering losses. We have thus assumed that the 
interaction of the mode fields with the fiber boundary causes signif- 
icant losses to high-order modes. On this basis we were able to calculate 
the fiber loss caused by random bends of the waveguide axis. For 
bends that approach a sinusoidal shape, with a period comparable to 

t We use (7i of (40) instead of <n of (42) since R was computed for q = 0. 
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the ray oscillation period in the parabolic index medium, the excess 
losses are extremely high. Bending of the waveguide axis with a period 
equal to the ray oscillation period must be avoided. For a statistical 
model, based on the assumption that the waveguide is composed of a 
sequence of circularly bent sections with random length and random 
radius of curvature, the waveguide losses have been predicted. We 
conclude that average radii of curvature of approximately 1 m can be 
allowed if an excess loss of 10 dB/km can be tolerated. The waveguide 
sections were assumed to have an average length of 1 cm. 

It is possible to reduce the width of the impulse response of a 
parabolic index fiber by coupling its modes by random bends of the 
fiber axis. The impulse response of parabolic index fibers is already 
quite favorable compared to the impulse response of the conventional 
optical fiber with a discontinuous but piecewise constant index dis- 
tribution. Our analysis shows that additional reduction of pulse delay 
distortion is accompanied by losses. A reduction of the pulse width 
to half its uncoupled width increases the loss by 0.7 dB, a ten-fold 
pulse width reduction increases the fiber loss by 18 dB. 
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APPENDIX 

The averaging process used to obtain (42) can be justified as follows. 
Each mode group characterized by the mode number q comprising 
all modes with 

< v < (J 2 n ° ka ~ 9 ) 

has the loss coefficient oi(q). By definition this can be written 

AP(q) is the power lost from the mode group per unit length and P(q) 
is the power carried by these modes. If we assume, for simplicity, that 
each mode carries the power P we can write P(q) = N(q)P so that we 
have 

ai{q) ~ WW 

with N(q) indicating the number of modes in the group. The total 
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loss is 



0"! = 



EAP(g) Y.<n(q)N(q) 

'i _ i 



E P(q) E N(q) 

'1 Q 

Replacing the sum by an integral yields formula (42). 
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